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ABSTRACT 

The fireball, the promising model of the gamma-ray burst (GRB), is an 
opaque radiation plasma, whose energy is significantly greater than its rest mass. 
We numerically simulate the evolution of the fireball heated by the neutrino- 
antineutrino annihilation process for the spherically symmetric case. We also 
derive analytical energy and momentum deposition rates via neutrino scatter- 
ing with thermalized electron-positron pairs in the fireball. In our simulation 
the matter is provided around the neutrinosphere before neutrinos start to be 
emitted, and the energy is injected during a finite period of time tdur- In the 
acceleration regime the matter shell is pushed from behind by radiation pressure. 
The Lorentz factor of the shell reaches the maximum value rj at r ~ T] 2 ct^ nT . 
After the fireball enters the coasting regime, the velocity distribution in the shell 
becomes very flat. The shell expansion rate dW/dr can be much smaller than 
n~ 2 . The runaway of temperature of the fireball due to neutrino scattering with 
electron-positron pairs does not occur in most cases. The energy deposition due 
to scattering is not so significant. 

Subject headings: gamma rays: bursts — hydrodynamics — neutrinos 



1. INTRODUCTION 



The relativistic fireball (Shemi and Piran 1990; Rees and Meszaros 1992; Meszaros and 
Rees 1993; Sari and Piran 1995; Sari et al. 1996) is one of the most promising models of 
gamma-ray bursts (GRBs). The sudden release of a large quantity of gamma-ray photons 
into a compact region produces a fireball that is an opaque radiation plasma whose energy 
is significantly greater than its rest mass. Piran et al. (1993) and Meszaros et al. (1993) 
investigated the hydrodynamics of fireballs and established basic behaviour in their evolution. 
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Most of the matter and energy is concentrated within a narrow shell. In the radiation- 
dominated phase, the rest mass and radiation energy density in the fireball behave as p oc r~ 3 
and e oc r -4 (r is the location of the shell), respectively. The bulk Lorentz factor 7 increases 
in proportion to r. After the fireball enters the matter-dominated regime, as long as the 
shell is optically thick to scattering, the densities behave as p oc r~ 2 and e oc r -8 / 3 , and 
the shell eventually coasts with 7 = 77 = E/Mc 2 , where M and E are the total baryon rest 
mass and the total energy of the fireball, respectively. A simple analytical estimate and 
hydrodynamical simulations indicate that the shell width expands as W ~ r/rf after the 
comoving energy density drops below the baryon rest mass density, while the shell width 
remains constant during the radiation-dominated phase. The actual behaviour of the fireball 
may be more complex depending on the energy injection process of the central engine. A 
detailed study of the fireball in various situations will be important to determine the character 
of the central engine in the future. 

The rapid temporal variabilities observed require that the GRB itself must arise from 
internal shocks within the flow. The shortest time scale is about 1 msec (Bhat et al. 1992). 
The internal shock is generated by the collision of two shells with different speeds arising 
from inhomogeneity of the fireball. In this paper we assume that the inhomogeneity in the 
fireball can be represented by multiple fireballs arising intermittently. 

Notwithstanding the very high energy phenomenon (~ 10 52 -10 54 ergs for spherical sym- 
metry), the baryon density in the fireball must be extremely small. In order to reproduce 
GRBs the final Lorentz factors of shells should be 100-1000. This is the famous baryon con- 
tamination problem and still remains unsolved. The central engine of GRBs is still beyond 
deep mist. The source of GRBs may be one super massive (failed) supernovae (Woosley 
1993; Paczyhski 1998) or may be a merger of two neutron stars or of a neutron star and a 
black hole (Eichler et al. 1989; Narayan et al. 1992; Meszaros and Rees 1992a; Katz 1997; 
Ruffert and Janka 1998, 1999). 

In these systems the hot core or accretion disk can emit neutrinos and antineutrinos. 
The neutrino-antineutrino annihilation into electrons and positrons (hereafter neutrino pair- 
annihilation) is a possible and important candidate to explain the energy source of GRBs 
(Paczyhski 1990; Meszaros and Rees 1992b; Janka and Ruffert 1996; Ruffert at al. 1997; 
Ruffert and Janka 1998, 1999; Salmonson et al. 2001). The central engine releases energy in 
a finite period of time at least longer than L/c where L is the scale of the central engine. 
Motivated by the delayed explosion of Type II supernovae, the energy deposition rate due 
to neutrino pair-annihilation above the neutrinosphere has been calculated (Goodman et al. 
1987; Cooperstein et al. 1987; Berezinsky and Prilutsky 1987). The energy deposition rate 
is proportional to r~ 8 for large r, and almost all deposition occurs near the neutrinosphere. 
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In hydrodynamic simulations of neutron star mergers, Ruffert and Janka (1999) showed that 
neutrino pair-annihilation deposits energy in the vicinity of the torus at a rate of (3-5) x 10 50 
ergs s -1 . In the failed supernovae model MacFadyen and Woosley (1999) predict that the 
total energy deposited along the rotational axes of Kerr black holes is (1-14) x 10 51 ergs. 

Most of the neutrinos' energy (> 99 %) will escape to infinity without annihilation. 
Since GRBs are very high energy phenomena, it is favorable that the energy of neutrinos 
is transformed into the fireball much more efficiently Asano and Fukuyama (2000, 2001) 
showed that the gravitational effect does not change the energy deposition rate dramatically. 
There is also another possible process to transfer the neutrinos' energy to the fireball. In the 
fireball, high-energy photons produce a large number of electron-positron pairs at the initial 
stage. The large optical depth due to electron-positron pairs leads to thermal equilibrium 
(Shemi and Piran 1990). The neutrino scattering with the thermalized electron-positron 
pairs (hereafter neutrino scattering) can inject energy into the fireball. This possibility was 
already pointed out by Woosley (1993). 

In this paper we discuss the evolution of the fireball heated by neutrinos and the for- 
mation of shells using a simple spherical model. We do not deal with a specific model of the 
central engine. We discuss only the general evolution of the fireball with energy injection 
in a finite period of time and the effect of neutrino scattering. While baryonic matter has 
been considered to be injected together with radiation energy in general, we assume here 
that the matter is provided around the central engine before neutrinos start to be emitted. 
In §2 we review the formulae for the energy and momentum deposition rates via neutrino 
pair-annihilation and derive the deposition rates via neutrino scattering. In §3 we numeri- 
cally simulate the fireball heated by neutrinos to investigate the evolution of the fireball and 
the deposition rates via neutrino scattering. Finally, §4 is appropriated for the summary. 

2. ENERGY AND MOMENTUM DEPOSITION RATE 
2.1. Pair Annihilation 



In this subsection we review the formulae for the energy and momentum deposition 
rates due to neutrino pair-annihilation. For simplicity, we assume that neutrinos are emitted 
isotropically from the neutrinosphere of radius R u . The neutrino number density in the phase 
space, f v d 3 p v = n u {e v )e 2 de u dVL v , is assumed to be conserved along the neutrinos' trajectory. 
This is a good approximation as long as the energy efficiency of neutrino pair-annihilation 
is very small. Defining X = [1 — (R^fr) 2 ] 1 / 2 , Goodman et al. (1987) derived the energy 
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deposition rate due to neutrino pair-annihilation as 
dE 



■ 12n 2 c 2 Rt 

pair v 



A|>G ' P N u N p ((e u )(el) + (el)(e D ))F E ( 



cdtdV 
where 

F s (r) = (l-X) 4 (X 2 + 4X + 5), (2) 

and the Fermi constant G| = 5.29 x 10 _44 cm 2 MeV~ 2 . The total emmision rate of neutrinos, 
N V) and the average energy of neutrinos, (e„), are written as 

N v (el) = (47r J R 2 )(7rc) | e^ k n v {e v )de v . (3) 

If we simplify the energy distribution of neutrinos by using the effective temperature, T e g, 
we can express the total emmision rate as 

where ((n) is the Riemann zeta function. The constant K p is written by the Weinberg 
angle sin 2 #vv = 0.23. To be exact we should sum up contributions of all species of neutrino. 
Each species may be emitted at different emission rates. When v c and v c are dominant 
components, K p = 0.124. If N u (e u ) and N v (e 2 u ) of all the six species of neutrino are the 
same, we can include the contributions of and v T by setting K p = 0.178. The fraction and 
mean energy of each species of neutrino depend on the model of the central engine (see, e.g., 
Ruffert at al. (1997)). In most cases the effective K p may be between 0.124 and 0.178 unless 
the emission rate of or u T is much larger than the rate of v c . For r ^> R u , the energy 
deposition rate behaves as oc T c 9 ff r~ 8 . Most of the energy is deposited near to the surface of 
the neutrinosphere. We also obtain the momentum deposition rate in a similar way as 



dP 



dtdV 

where 



A N,N,((e,)(el) + (el)(e,))F P (r), (5) 



F P {r) = (1 - X) 4 (l + X)(3X 2 + 9X + 8)/4. (6) 



2.2. Neutrino Scattering 



In the fireball high-energy photons produce a large number of electron-positron pairs. 
The large optical depth due to electron-positron pairs leads to thermal equilibrium (Shemi 
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and Piran 1990). Neutrino scattering with thermalized electron-positron pairs can transfer 
the neutrinos' energy to the fireball. We derive the energy and momentum deposition rates 
through neutrino scattering in this subsection. The electron-positron pairs are assumed to 
be relativistic e c ^> m c c 2 . For the reaction u e + e _ — > v c + e~, the differential cross-section 
is given by 

^ = ^ fa 2 ; P *' c) [(2 sin 2 9 W + I) 2 + 4 sin 4 9 W (1 - yf] , (7) 

where y = p ijB ■ (p ijV - Pf, e )/(Pi,v • Pi,e) is Lorentz invariant, and (pi, u -Pi, e ) expresses the inner 
product of the 4-momenta of a neutrino and an electron. The subscripts "i" and "/" denote 
the initial and final states, respectively. We can write the cross-sections for the other five 
species of neutrino in a similar way (see, e.g., Kim and Pevsner (1993)). 

In the comoving frame of the relativistic fireball, the energy of the neutrinos redshifts 
as e' v = £„/7(l + /3cos0J,), where 7 = 1/(1 — (i 2 ) 1 / 2 is the Lorentz factor of the fireball and 9' v 
is the incident angle of neutrinos to the r— axis (see Fig. 1). Hereafter values with a prime 
are measured in the comoving frame. We have two scattering angles for a given value of y. 
The energy a neutrino loses per one scattering depends on the final scattering angle. The 
average energy an electron/positron obtains for a given value of y becomes As' = y{s' u — e' e ). 
The average momentum component an electron/positron obtains, which is parallel to the 
neutrino's momentum, becomes cApj| = y{e' v — e' c cos9'), where 9' is the incidence angle 
between the two particles. Here we assume that the random motion of electron-positron 
pairs is isotropic in the fireball frame. When a neutrino interacts with electrons distributed 
isotropically, the mean vertical momentum the neutrino loses per one scattering is zero. As 
we have assumed in §2.1, if n v (e v ) of all the six species of neutrino are the same, the energy 
a scattered electron (positron) gains and the total cross section are averaged as 



n us (e' u )a s ■ As' = n u (e' u ) ^ ~^ y ^ ~ £ '^ dy 



(8) 



= ^n^)K s G 2 M,v-p'i, e W v -<), (9) 

where the summation is over all species of neutrino, and K s = (24 sin 4 6* w — 4 sin 2 9 w +3)/2n = 
0.533. When v c and P e are dominant components, K s should be altered as (16 sin 4 #w + 
Ssin^w + 2)/47r = 0.373. To be more precise we have to estimate the deposition rate 
corresponding to the variation of the distributions of the six species of neutrino. We include 
the variation by adjusting K s effectively as we have done in §2.1. 



The numbers of neutrinos and antineutrinos are assumed to be the same. Then the 
reaction rates of electrons and positrons are the same because of the symmetric property. 
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Since the energy distribution of the positrons is the same as the distribution of the electrons 
in thermal equilibrium, we include the contribution of positrons by doubling the energy 
distribution of electrons n e (e' e ). The energy deposition rate in the comoving frame is written 

as 



dE 



cdtdV 



= I I I /^<^Q^£%V S (4)^e(<)^-A £ 'c(l-COS^)(lO) 
scat J J J J ^ 

J J de' u de' e n„(e'„)n e (e' e )e*ef(e' u - e' e ) 

t / d0'sin0'(l -cos9') 2 , (11) 
Jo 



UcK s G 2 F ir 2 



J 

Jo 



max 

/ d0' u sin 9'„ 



where cos6^ ax = (X — — f3X), including the redshift effect. Assuming the electron- 

positron pairs are in thermal equilibrium and their temperature T c ^> m c c 2 , we get 



dE ' 
cdtdV 



f(T c (r))G' E (T e ,(3;r), 



scat 



where 



(12) 



(13) 



and 



G> E (T e ,P;r) = 7 2 (1-X){ 7 ^(e 2 ) 



1--/3(1 + X )+ 
(3 2 (1 + X + X 2 )- ^-(1 + X)(l + X 2 ) 



-20C(5)(e,)T e (r) 



1-/3(1 + X) + ^(1+X + X 2 ) 



(14) 



The temperature T e (r) is measured in the comoving frame. The momentum deposition rate 
can be obtained in a similar fashion, using the parallel component of momentum to the 
r-axis, cos 9' u Apj| . Carrying out the Lorentz transformation, we obtain the energy deposition 
rate in the observer frame as dE/cdtdV\ scait = /(T e (r))G s (T e , (3; r), where 



G £ (T c ,/3;r)= 7 { 7 ^( e 2 > 
-20C(5)(e,)T e (r) 



l-[3(l + X) + ^(l + X + X 2 ) 
i(3 7 2 -l) + ^(l- 6 7 2 )(l+X) 



+\{i 2 - 1)(1 + X + X 2 ) 



[l-X). 



(15) 
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In the case T e (r) < m e c 2 these formulae overestimate the rates. These equations give the 
upper limit of the energy and momentum deposition rates. For (3 ~ and r ^> R u , in- 
troducing the effective temperature of the neutrinos, the energy deposition rate behaves as 
oc (T c 5 ff T c 4 (r) — (r)T^ s )r~ 2 . The deposition rate depends on the temperature of the fire- 
ball. When r = R v and f3 — 0, the constant factor, except for the temperature dependence, 
in the deposition rate via neutrino scattering is larger than the factor via pair-annihilation 
by a factor of ~ 5. If T e ~ T eff , the energy injection rate for neutrino scattering becomes 
comparable to the rate for neutrino pair-annihilation. However, if T c > T eS the scattering 
effect cools down the fireball inversely. The deposition rate via pair-annihilation behaves as 
oc r~ 8 . If T c 4 (r) has a flatter distribution than r -6 , much more energy could be injected by 
scattering than energy injected by annihilation at a great distance. 

The momentum deposition rate in the observer frame is obtained by dP/dtdV\ scat = 
f(T e (r))G P (T e ,(3;r), where 



G P (T c ^-r)= 1 { 1 ^{el) 



' " ' ' 01 1 '1 + X) - 1(3(1 + X + X 2 ) + ^(1 + X)(l + X 2 ) 



+ 10((5)(e u )T e (r) 



2 V ' 3 v 4 
v _ , X)-3/3 7 2 -^(l + 2/3 2 )(l + X + X 2 )j|(l-X).(16) 



For f3 ~ 0, neutrinos accelerate the fireball regardless of T e . Neutrino scattering can decel- 
erate the flow for large 7 like Compton dragging. 

There is a possibility that the temperature of the fireball could run away via neutrino 
scattering (Woosley 1993) and reach T e ~ T c g. In addition the total energy could be enhanced 
by the energy deposition at large r if T e (r) distributes flatter than r~ 6 . However, it seems 
difficult to lead T c to the runaway, because the dependence ~ T e 4 makes the deposition 
rate due to neutrino scattering much weaker than the rate due to neutrino pair-annihilation 
giving a slightly smaller T e than T eff . To reach higher temperatures the expansion speed of 
the plasma should be slow enough during the energy injection such that the heating effect 
via neutrinos is dominant rather than adiabatic cooling. In most cases we need to simulate 
the evolution of the fireball numerically in order to estimate the conclusive energy deposition 
via neutrinos. 



3. NUMERICAL CALCULATION 

In this section we numerically simulate the fireball heated by neutrinos. Neutrinos are 
assumed to be emitted isotropically from the neutrinosphere during a finite period of time. 
As Piran et al. (1993) assumed in their calculation, we assume the fireball is optically thick. 
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The radiation and matter at each radius behave like a single fluid moving with the same 
velocity. Since the radiation pressure dominates, the pressure p and the energy density e 
of radiation and relativistic electron-positron pairs are related by p = e/3. The relativistic 
conservation equations of baryon number, energy and momentum are written by 



ld_ 

cdt 



(pi) 



]_d_ 

r 2 dr 



(r 2 pu) , 



(17) 



1 d /4 



cdt U 67 



]_d_ 

r 2 dr 



r ' 2 l u | g e + P 



+ 



dE 
cdtdV 



(18) 



1 d 




1 d 




cdt 


U1 + P ) 




r 2 u 2 + p^j 



lde dP 
+ 



3dr dtdV 



(19) 



where u = 7/3. The rest mass energy density of baryons p and the non-baryonic energy 
density e are measured in the comoving frame. 

We use the relativistic "Roe-solver" method in the computation of one-dimensional rel- 
ativistic flow. The Riemann problem is approximated numerically by means of a specially 
chosen linearized form, requiring less computational effort. This Eulerian code was devel- 
oped for the relativistic case by Mellema et al. (1991). However, our computation is not 
so complicated as the computation of Mellema et al. (1991), because we assume spherical 
symmetry and neglect the gravitational effect. Our code is relatively simple and has passed 
several standard tests. In addition we have checked our numerical code by reproducing the 
results of the simulation in Piran et al. (1993). We adopt the spatial and time resolution 
as Ar = i^/3000 and cAt = Ar/10, respectively. In order to discuss the evolution of 
shell width, however, we need to simulate with a very high spatial resolution. The expan- 
sion/contraction rate of shell width may be limited as \dW/dr\ < I/7 2 . It will take very long 
time to follow the change of shell width accurately for large 7 because of the high resolution. 
In our simulation, judging from test calculations, the expansion/contraction rate is not so 
reliable for 7 > 10. However the rate of change is small enough for large 7, in our calculation, 
such that we can follow the approximate behaviour of the fireball. 

First we investigate the evolution of the fireball in an ideal situation to produce a 
relativistic flow. Cold baryonic matter is provided around the neutrinosphere before the 
simulation starts. Then we study the effect of neutrino scattering on the energy deposition 
rate. 
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3.1. Evolution of the Fireball 

In order to produce a relativistic flow, we assume an ideal situation. The initial profile 
of the matter that we present is compact as p(r,t = 0) oc (Rq + r 8 ) -1 for r > R v . The 
matter is initially cold and static as e(r, 0) = and 7(7", 0) = 1. In this simulation the energy 
is injected at the constant rate, E, during < t < (r — R u )/c + t dur , while in previous 
studies they had provided radiation energy before the simulations started (Piran et al. 1993; 
Meszaros et al. 1993). Spatial distribution of E is given by equation (2). We neglect the 
effect of neutrino scattering in this subsection. We set the initial core radius of the matter 
as i? = ctdur hereafter. The rest mass energy is normalized as 4nr] J drr 2 p = Et^nr, where 
77 is the ratio of the radiation energy to the rest mass energy. The mean final Lorentz factor 
of the matter shell should be 77 + 1. Thus, the physical conditions in this simulation are 
determined by the two parameters i] and Td ur = cta m /R u . 

First of all we show results for the case of rj = 60 and 

^"dur — ^ &n example. If we 
introduce physical units as R u = 10 km and T e g- = 10 MeV, the initial matter density at 
r = R u corresponds to 2.4 x 10 4 g cm -3 in this case. Figure 2-4 are the profiles of the 
fireball at different times. At r = ct/R u = 3 (Fig. 2), namely, at the end of the energy 
injection, the matter is swept out by the radiation pressure. The information about the 
initial configuration of the matter is almost lost. The approximation p < e, which has been 
assumed conventionally in discussions about the initial evolution of the fireball, already 
failed around the matter shell at this time. The matter shell is compressed by the radiation 
pressure. The FWHM of the matter shell in the observer frame (the profile of pj) is about 
0.1R U and remains nearly constant within the numerical error limit hereafter. The behavior 
of the fluid around the neutrinosphere agrees with the stationary solution of pure radiation 
flow with energy injection. We have confirmed that e and 7 behind the matter shell in Figure 
2 are very close to the numerical solutions obtained by setting the left-hand side of equations 
(17)-(19) to zero and p = 0. Figure 3 is the fireball at r = 6. The energy injection has already 
finished. The pure radiation flow, which had been deposited later, is chasing the matter shell. 
Around the matter shell the forward and reverse shocks are formed. The reverse shock is 
radiation dominant, while the forward one is matter dominant. In Figure 4 we display the 
fireball at r = 30. The Lorentz factor becomes a maximum behind the shell, where the 
radiation is dominant. The radiation flow pushes the matter from behind. The matter shell 
itself accelerates more slowly than the prediction by the simple analysis, 7 = (r/R u ). The 
behaviour of u at the point where p7 has the maximum value is approximated as um ~ ^°' 5 
at this moment. The value um does not accelerate obeying a simple power law. The power 
law index declines as um oc r 4 , r°- 3 ,... with increasing r. We need a very long time to reach 
Mm = i] in our computation. The maximum values of p and e do not obey simple power-law 
behaviour. 
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In Figure 5 we plot um for different values of Td ur and r] = 60. The acceleration time 
is roughly in proportion to the duration time of the energy release. When the shell reaches 
r = 1]R U , 7 becomes a few tens of a percent of 77. The maximum Lorentz factor in the whole 
system increases in proportion to r. Extrapolating um from Figure 5, we guess um may reach 
i] at r ~ i] 2 TdurRu- At the end of the energy injection the distance from R v to the head of 
the shell is about t^Rv Therefore, the energy injected at the end will catch up with the 
matter shell of 7 ~ rj at r ~ rfr Am R u . This is the reason u M reaches 77 at r ~ rfT Aur R v . 

To economize the calculation time we demonstrate the evolution of the fireball for i] = 10 
and Tdur = 1. In Figure 6 we plot um and FWHM of the profile of fry. Since the initial matter 
distribution is narrow in comparison with the case r dur = 3, the radiation pressure expands 
the matter component initially. Then the matter component is compressed. Figure 7 shows 
that the rear part of the matter shell is faster than the head part at r = 12 that implies the 
shell is compressed. The shell begins to expand again from r ~ 30. In this period the rear 
part of the matter shell is slightly slower than the head part, as is shown in Figure 8. The 
expansion rate is dW/dr ~ 10 -3 < rf 2 , which is close to the spatial resolution Ar/R u in 
our simulation. This value may not be accurate because of the size of the numerical error. 
However, the difference of the Lorentz factor in the shell is rather small as £7 <C 7. At 
t = 105 the maximum and minimum u within FWHM are 10.28 and 9.89, respectively. The 
maximum relative velocity within FWHM in the observer frame is 3.8 x 10~ 4 c that suggests 
dW/dr ~ 4 x 10~ 4 , which is smaller than our error limit in the computation. Thus, we can 
expect a smaller expansion rate than rf~ 2 by a factor of 10 or 10 2 . As we have assumed, um 
reaches rj at r/R u ~ r ~ if = 100. We can expect that the behaviour of the fireball for large 
77 would be similar to the above results. 

We refer the simulation in Piran et al. (1993) to compare with the case in which radiation 
energy has been given before the simulation starts. In their simulation the initial profiles 
are given as e(r, 0) oc (R* + r 8 ) -1 and p(r, 0) = e/rj. According to Figure 2 in Piran et al. 
(1993) the matter shell expands even in the acceleration regime. At r ~ 100 in our notation 
the shell width grows to a few times the size of Rq. The Lorentz factor in the matter shell 
increasingly distributes with radius and their difference in the shell is rather large as £7 ~ 7. 

The results of this subsection can be summarized as follows: In the acceleration regime 
the shell is pushed from behind by the radiation pressure. The expansion of the shell is 
suppressed in this regime. The acceleration declines as the shell evolves. The Lorentz factor 
of the shell reaches 77 at r ~ i] 2 Td UT R u . After the fireball enters the coasting regime, the 
shell begins to expand. The velocity distribution in the shell is very flat. We need a more 
accurate numerical simulation and longer calculation time to estimate the expansion rate. 
There is a possibility that the expansion rate is much less than r\~ 2 . The simulation in 
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this subsection is merely one example of the evolution of the fireball being affected by the 
energy deposition due to neutrinos. However, when radiation fluid is deposited inside the 
initial matter distribution during a finite period of time, the fireball may evolve in a similar 
manner regardless of the details of the spatial distribution of E and initial matter. 

3.2. Effect of Neutrino Scattering 

In this subsection we simulate the fireball including the effect of neutrino scattering. 
Our purpose is to investigate the enhancement of the energy deposition by the scattering 
effect. The energy deposition by neutrinos including the scattering effect is a nonlinear 
process. The rate of neutrino scattering depends on the temperature of the fireball. As is 
shown in the former subsection, the temperature around the neutrinosphere can be expressed 
by the stationary solution of the pure radiation flow as long as the initial matter density is 
sufficiently low. It is easily found that the maximum temperature, which is obtained from 
the stationary solution, is much lower than the effective temperature of the neutrino. In an 
ideal situation for producing the relativistic shell as in §3.1, the effect of neutrino scattering 
is negligible. 

However the environment is generally very dense around candidates of the central engine: 
the central region of massive stars, merging neutron stars, and so on. In a very dense 
environment the matter cannot be swept out promptly. The energy injected by neutrinos 
may stand near to the neutrinosphere, and the fireball temperature becomes higher than the 
case we considered in the former subsection. We investigate the effect of neutrino scattering 
for a dense environment. In order to consider the most effective case on neutrino scattering, 
we make the matter density constant as 10 9 g/cm 3 and set R u = 10 km and T e g = 15 
MeV. Similarly, as in the former subsection, the energy deposition due to pair-annihilation 
is given by equation (1) with K p = 0.178. We add the contribution of the scattering effect 
discussed in §2.2 to equations (17)-(19) with K s = 0.533. The temperature of the fireball 
is obtained from T e 4 = 60(ftc) 3 e/(ll7r 2 ). Since our purpose is to investigate the effect of 
scattering between neutrinos and electron-positron pairs in the extreme environment, we 
ignore the interaction between baryons and neutrinos. As a matter of fact, for T c g = 15 
MeV, the energy deposition rate due to the neutrino-baryon interaction (oc T c 6 ff p for f3 — 0) 
becomes dominant rather than the rate due to pair-annihilation at r = R u , in a case in 
which the initial matter density is larger than 7 x 10 9 g/cm 3 . As long as p is constant, the 
total energy deposition via the neutrino-baryon interaction is larger than the deposition via 
pair-annihilation. 

In Figure 9, we display the fireball at r = 1. We have normalized the length scale and 
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energy density by 10 km and (10 5 g cm~ 3 )c 2 , respectively. As we have assumed above, the 
expansion speed is slow and the radiation energy stands inside the shell. The distribution of 
e inside the matter shell is nearly flat. The typical temperature T e inside the matter shell, 
where the radiation is dominant, is about 5 MeV, smaller than the effective temperature of 
neutrinos. At r = R v the scattering effect enhances the deposition rate by only 5 %, while 
the total enhancement is 20 % at r = 10 (see Fig. 10). This implies that a considerable 
amount of energy is deposited in a wide range of r by scattering. 

The temperature T C (R V ) increases after the start of the energy deposition (see Fig. 11). 
However, the huge energy injection (T e g = 15 MeV) sweeps out the matter immediately. 
Then adiabatic cooling due to the expansion becomes stronger than the neutrino heating. 
The temperature T e around the neutrinosphere declines and reaches the value given by the 
stationary solution. 

If we enlarge the dimension of the region as R v = 100 km, the energy deposition rate 
(oc Rl) and the energy efficiency (oc R v ) due to pair-annihilation increase. The dynamical 
time scale becomes larger by a factor of 10, while the energy injection rate per volume 
remains constant for the given value of r/R u . Therefore, the injected energy per baryon 
during a dynamical time scale increases, which leads to high temperatures. In our case the 
fireball temperature T e increases to ~ 7 MeV (Fig. 11). The energy deposition is enhanced 
by 60% because of scattering at r = 10. The qualitative behaviour in this case is almost the 
same as in the former case. 

On the other hand, if the effective temperature decreases to T cS = 5 MeV, the smaller 
energy injection cannot promptly sweep out the matter. The adiabatic cooling effect is 
weaker than the heating effect, initially. The temperature T C (R U ) increases with time, but 
remains much smaller than T e g- (Fig. 11) in comparison with the case T c g = 15 MeV. Then 
the increasing rate of T C (R V ) is saturated. To make matters worse, the energy deposition 
due to the neutrino-baryon interaction is practically dominant for lower T e . Neutrino pair- 
annihilation and scattering play secondary roles. Such a situation is not adequate for our 
purpose. Reduction of the initial matter density leads to a prompt sweeping out of matter. 
Thus, it is very difficult to make T e large. 

The temperature of the fireball does not run away in our cases and the nonlinear effect 
is not so strong. The energy deposition due to scattering is not so important even under the 
extreme condition we adopted. 
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4. SUMMARY 

We have simulated the fireball in which radiation energy is deposited in a narrow region 
(E ~ r -8 ) during a finite period of time. The baryonic matter is initially prepared around 
the neutrinosphere. In the acceleration regime the matter shell is pushed from behind by the 
radiation pressure. The Lorentz factor of the shell reaches rj at r ~ rj 2 T dui R u . Meszaros and 
Rees (2000) estimated that wind becomes optically thin at ~ 10 6 i?^ with typical parameters. 
If rd ur 3> 1, the shell can be optically thin before 7 reaches r\ even if 77 is smaller than the 
critical value 77* (see Meszaros and Rees (2000)). The maximum values of p and e do not 
obey a simple power-law behavior. After the fireball enters the coasting regime, the shell 
may begin to expand. The velocity distribution in the shell is rather flat in comparison with 
the case in which the radiation energy has been given before the simulation starts. Although 
we need a more accurate numerical simulation to estimate the expansion rate quantitatively, 
dW/dr may be smaller than rf 2 . 

The short time scale of the temporal profile of GRBs requires that the shell width 
should remain thin at the collision radius. If the shell width expands as W = r/rj 2 , the 
width becomes ~ 10 10 cm for rj = 100 at the typical collision radius ~ 10 14 cm. Then 
the observed emission time scale of one collision becomes ~ 1 sec, which is longer than the 
observed timescale of the pulse rise. In addition, if the shell expands, the magnetic field in 
the shell may weaken. It is very difficult to reproduce hard spectra of observed GRBs in this 
case (Guetta et al. 2001). Therefore, it is favorable that the shell width is constant. The 
actual shell expansion depends on the Lorentz factor distribution across the shell, which is 
determined by the energy ejection process around the central engine. Although the actual 
activity of the central engine is more complex than our case, our simulation shows a possibility 
that the expansion rate dW/dr is smaller than i]~ 2 . Our simulation is one simple example 
of the evolution of the fireball. If baryonic matter is also injected with radiation energy, the 
behaviour of the fireball may be close to the behaviour in Piran et al. (1993). In our case 
baryonic matter should be injected during the quiescent period of the central engine. The 
matter could be injected from the side if the central engine has a jetlike structure. But we 
do not know if this is realistic or not. In the future, we will have to examine whether the 
shell width remains narrow in some kind of model of the central engine. Future observations 
may constrain shell's structure and give a clue to determining the central engine. 

As for the effect of neutrino scattering, the temperature of the fireball does not run 
away in most cases that we assume, and the nonlinear effect is not so strong. Although 
the scattering effect can enhance the energy deposition rate by a few tens of a percent, the 
effect on the deposition rate is not so important even for the extreme conditions we assumed. 
When other processes of energy injection, like magnetic reconnection, etc., provide a large 



-14- 



number of pairs, the effect of neutrino scattering can be larger. In any case a very dense 
environment is required for raising the fireball temperature high enough in a wide range 
of r. Under the dense environment the baryon contamination is crucial, while the higher 
temperature increases the energy deposition due to the scattering effect. Since the energy 
enhancement due to the scattering effect is not significant, it is difficult to make 7 large 
enough in the situation where the effect of neutrino scattering becomes important. 

We thank Dr. A. Iida for his useful advice concerning the hydrodynamical simulation. 
The authors are supported by a Research Fellowship of the Japan Society for the Promotion 
of Science. 
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Fig. 1. — Definitions of 9 max , V , and 6. 
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Fig. 2. — Fireball at r = 3 for r] = 60 and Td ur = 3. The thick and thin solid lines are fry 
and ej, respectively. The dotted line is u = 7/3; p7 and are plotted in logarithmic scale. 
Here we neglect the effect of neutrino scattering. 
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Fig. 5. — Values of %: the value u at the maximum point of p7 for 
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Fig. 8. — Same as Fig. 7, but r = 105 
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Fig. 9. — Fireball at r = 1 in the very dense case (see §3.2). The initial matter density is 10 £ 
g cm -3 , T cff = 15 MeV, and R u = 10 km. Here we include the effect of neutrino scattering. 
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Fig. 10. — Total energy deposited by neutrinos (solid line), including the effect of neutrino 
scattering and the energy deposited via neutrino pair-annihilation only (dotted line). The 
physical parameters are the same as those in Fig. 9. 
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Fig. 11. — Fireball temperature T c at r = R v . The initial matter density is 10 9 g cm~ 3 in 
all cases. The effective temperature T e g and R v are indicated. 



